%nonsingular part of core

function y = func_NSpart(zeta, m, CONSTS)

    Sm     = func_Sm(zeta, m, CONSTS);
%     if(m==0)
%         Sm_mod = func_Sm_mod(zeta, m, CONSTS);
%     elseif(m==-1)
%         Sm_mod = func_Sm_mod(zeta, m, CONSTS);
%     else
%         Sm_mod = 0;
%     end
%     if(m==-1)
        Sm_mod = func_Sm_mod(zeta, m, CONSTS);
%     else
%         Sm_mod = 0;
%     end

    y = Sm + Sm_mod;

end